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ABSTRACT 

We present a deep X-ray observation of the low-luminosity active galactic nucleus in NGC 4258 (M 106) using 
the Advanced Satellite for Cosmology and Astrophysics (ASCA). Confirming previous results, we find that the 
X-ray spectrum of this source possesses several components. The soft X-ray spectrum (< 2keV) is dominated 
by thermal emission from optically-thin plasma with kT 0.5 keV. The hard X-ray emission is clearly due to 
a power-law component with photon index F w 1.8 absorbed by a column density of Ah ~ 8 x 10^^ cm^'^ . 
The power-law is readily identified with primary X-ray emission from the AGN central engine. Underlying both 
of these spectral components is an additional continuum which is possibly due to thermal emission of a very hot 
gaseous component in the anomalous arms and/or the integrated hard emission of X-ray binaries in the host galaxy. 
We also clearly detect a narrow iron Ka emission line at ^ 6.4 keV. No broad component is detected. We suggest 
that the bulk of this narrow line comes from the accretion disk and, furthermore, that the power-law X-ray source 
which excites this line emission (which is typically identified with a disk corona) must be at least ~ 100 GM/c^ 
in extent. This is in stark contrast to many higher-luminosity Seyfert galaxies which display a broad iron line 
indicating a small (^10 GM/c^) X-ray emitting region. It must be stressed that this study constrains the size of 
the X-ray emitting corona rather than the presence/absence of a radiatively efficient accretion disk in the innermost 
regions. If, instead, a substantial fraction of the observed narrow line originates from material not associated with 
the accretion disk, limits can be placed on the parameter space of possible allowed relativistically broad iron lines. 
We include a discussion of various aspects of iron line limb-darkening for highly inclined sources, including the 
effect of gravitational light-bending on the apparent limb-darkening law. By comparing our data with previous 
ASCA observations, we find marginal evidence for a change in absorbing column density through to the central 
engine, and good evidence for a change in the AGN flux. We conclude with a brief discussion of two serendipitous 
sources in our field of view; the QSO Q1218H-472 and a putative z ~ 0.3 cluster of galaxies. 

Subject headings: accretion, accretion disks - black hole physics - galaxies: active - galaxies: individual: 
NGC 4258 - X-rays:galaxies 



1. INTRODUCTION 

The nearby low-luminosity active galactic nucleus (AGN) in 
NGC 4258 (M 106) has become crucially important in our un- 
derstanding of accreting supermassive black holes. Exquisite 
position and velocity measurements of the H2O megamasers 
in NGC 4258 reveal that the masing material resides in a very 
thin and slightly warped disk that is in almost perfect Keplerian 
motion about a central black hole with a mass of 3.5 X 10"^ Mq 
(Nakai, Inoue & Miyoshi 1993; Greenhill et al. 1995a; Miyoshi 
et al. 1995). Detailed studies of the maser proper motions and 
centripetal accelerations confirm this interpretation (Greenhill 
et al. 1995b) and allow us to measure the distance to this source 
independently of the Hubble constant (d = 7.2 Mpc; Herrn- 
stein et al. 1999). 

Sensitive X-ray observations provide a powerful means of 
probing both large scale and small scale structures within 
NGC 4258. X-ray emission was first detected in a short Ein- 
stein observatory high resolution imager (HRI) observation 
(Fabbiano et al. 1992). More sensitive soft X-ray data col- 
lected by the ROSAT position sensitive proportional counters 
(PSPC) and HRI found an extended halo of hot (4 x 10^ K) gas 
around NGC 4258 (Pietsch et al. 1994; Vogler & Pietch 1999) 
as well as X-ray emission associated with the well known he- 
lically twisted jets (the anomalous arms; Pietsch et al. 1994; 
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Cecil, Wilson & De Pree 1995; Vogler & Pietch 1999). None 
of these soft X-ray observations penetrated the column of ab- 
sorbing gas that obscures the AGN itself. This had to await 
ASCA observations (Makishima et al. 1994; hereafter M94). 
M94 found that the central X-ray source was well described 
by a power-law with photon index F « 1.8, and was absorbed 
by a column density of Ah 1.5 x 10^^ cm^^ . The soft X- 
ray spectrum was found to be complex with components aris- 
ing from thermal plasma emission and, possibly, contributions 
from the underlying X-ray binary population in NGC 4258. A 
marginal detection of an iron line was claimed with an equiva- 
lent widthof 250± lOOcV. 

Of great importance is the opportunity that NGC 4258 gives 
us to study accretion physics when the mass accretion rate is po- 
tentially very low. The significance of this issue is highlighted 
when it is realized that most of the (quiescent) supermassive 
black holes in the universe are inferred to accrete matter at rates 
which are comparable to, or less than, that in NGC 4258. It was 
realized by several authors that when the accretion rate is low 
(relative to the Eddington rate), an accretion disk may switch 
into a hot, radiatively-inefficient mode (Ichimaru 1977; Rees 
1982; Narayan & Yi 1994; Narayan & Yi 1995). In essence, the 
plasma becomes so tenuous that the timescale for energy trans- 
fer from the protons to the electrons (via Coulomb interactions) 
becomes longer than the inflow timescale. The energy remains 
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as thermal energy in the protons (which are very poor radia- 
tors) and gets advected through the event horizon of the black 
hole. These are the so-called Advection Dominated Accretion 
Flows (ADAFs). ADAFs are to be contrasted with 'standard' 
radiatively-efficient accretion disks in which the disk remains 
cool and geometrically thin all of the way down to the black 
hole (Shakura & Sunyaev 1973; Novikov & Thorne 1974). 

X-ray observations of broad iron Ka lines have shown that 
higher luminosity systems do indeed accrete in the radiatively- 
efficient mode (Tanaka et al. 1995; Fabian et al. 1995). How- 
ever, the basic nature of the accretion disk when the mass ac- 
cretion rate is low is still far from clear. The existence of the 
ADAF solution is at the mercy of poorly known physics such 
as the strength of the electron-ion coupling and the fraction of 
the viscous energy that is deposited into the electrons (Quataert 
& Gruzinov 1999). Also, Blandford & Begeknan (1999) have 
suggested that the ADAF solutions discussed by Narayan & Yi 
(1994, 1995) are physically inconsistent and necessarily drive 
powerful outflows (producing the so-called Adiabatic Inflow- 
Outflow Solutions; ADIOS). Even if ADAF-type solutions ex- 
ist, it remains an open question as to whether real disks can 
make a transition to this mode when the outer regions of the 
disk are cold and geometrically-thin (as they clearly are in sys- 
tems such as NGC 4258). NGC 4258 provides a laboratory in 
which we can examine all of these issues. 

This dichotomy in possible accretion disk physics has pro- 
duced two models for the central regions of NGC 4258. By 
modeling the observed water maser emission in this system, 
Neufeld & Maloney (1995) conclude that the accretion disk in 
NGC 4258 has a high efficiency (;^ 10%), and a low accretion 

rate {M /a ^ 10^'* M0 yr""'^ , where a is the standard viscosity 
parameter of accretion disk theory). Furthermore, the observed 
maser emission in part traces out a warp in the disk. Model- 
ing this warp as being driven by radiation pressure from the 
central X-ray source (Pringle 1996) also implies a radiative ef- 
ficiency ^ 10% (Maloney, Begelman, & Pringle 1996). On the 
other hand, Lasota et al. (1996) use ADAF models of the 2-10 
keV X-ray power law slope and continuum radio data to argue 
that the system has low efficiency (^ 0.1%) and high accretion 
rate {M/a - lO^^Mgyr^^ a ^ 0.3). While the original 
model of Lasota et al. (1996) postulated a large ADAF region 
(r ~ 10^ GM/c^), recent radio data constrain the ADAF region 
to be smaUer than r ~ 100 GMjc^ (Herrnstein et al. 1998). 

In this paper, we report a deep (^ 200ksec of good expo- 
sure time) X-ray observation of NGC 4258 with ASCA (Ohasi 
et al. 1996; Makishima et al. 1996; Yamashita 1997). We also 
obtained simultaneous hard X-ray data with the RXTE satellite 
but, due to a gain change prior to our observation, there are cur- 
rently no robust response matrices or background models for 
these particular RXTE data. A presentation of these RXTE data 
must await these developments in calibration. In section 2, we 
describe the data reduction before discussing the X-ray spec- 
trum of NGC 4258 in Section 3. Section 4 focuses on the prop- 
erties of the observed iron line and the implications of this line 
for the nature of the accretion disk. Section 5 presents a brief 
discussion of the other interesting objects in the field of view 
of the ASCA Gas Imaging Spectrometer (GIS). Our results are 
summarized in Section 6. 

2. DATA REDUCTION AND THE X-RAY IMAGE 

The NGC 4258 field was observed by ASCA on 1999 May 
15-20. The SIS data were cleaned in order to remove the ef- 



fects of hot and flickering pixels and subjected to the follow- 
ing data-selection criteria : i) the satellite should not be in the 
South Atlantic Anomaly (SAA), ii) the object should be at least 
7 degrees above the Earth's Umb, iii) the object should be at 
least 25 degrees above the day-time Earth limb and iv) the lo- 
cal geomagnetic cut-off rigidity (COR) should be greater than 
6 GeV/c. We also apphed a standard grade selection on SIS 
events in order to further minimize particle background. The 
GIS data were cleaned to remove the particle background and 
subjected to the following data-selection criteria : i) the satellite 
should not be in the SAA, ii) the object should be at least 7 de- 
grees above the Earth's limb and iii) the COR should be greater 
than 7 GeV/c. SIS and GIS data that satisfy these criteria shall 
be referred to as 'good' data. 

After the above data selection, there are 170ksec of good 
data per SIS detector and 185 ksec of good data per GIS detec- 
tor. Images were then extracted from these good data for each 
of the four instruments (two SIS and two GIS). 

Fig. 1 shows the GIS2 image for this observation. The nu- 
cleus of NGC 4258 is the brightest X-ray source in this field. 
Two other sources are also detected: (a) one point-like source 
~ 7 arcmins to the west of NGC 4258, and (b) another slightly 
extended source ^ 17 arcmins to the north-west of NGC 4258. 
Both of these sources were clearly detected and studied with 
ROSATby Pietsch et al. (1994). Source (a) can be readily iden- 
tified as file ^ « 0.4 quasar Q1218H-472 (Burbidge 1995; Bur- 
bidge & Burbidge 1997). The identification of source (b) is less 
secure, but Pietsch et al. argue that it is a background (z ^ 0.2) 
cluster of galaxies on the basis of a possible galaxy over-density 
on a deep optical plate. 

We have extracted spectra and lightcurves for aU three 
sources in our field. Unless otherwise stated, source counts 
were extracted from a circular region centered on the source 
with radii of 3 arcmin and 4 arcmin for the SIS and GIS respec- 
tively. Background spectra were obtained from blank regions 
of the same field (using the same chip in the case of the SIS). 
No temporal variabiUty was observed in any of these sources. 
In order to facilitate spectral fitting, all spectra were re- 
binned so as to contain at least 20 photons per spectral bin. 
In order to avoid poorly calibrated regions of the spectrum, the 
energy ranges considered were 1.0-lOkeV for the SIS detec- 
tors, and 0.8-lOkeV for the GIS detectors. Note that we use a 
lower-energy cutoff for the SIS that is considerably higher than 
the 'standard' 0.6 keV cutoff in order to avoid the effects of 
"residual dark current distribution", or RDD, which is known 
to plague recent ASCA observations. 

3. A DETAILED X-RAY STUDY OF NGC 4258 

We now discuss the X-ray properties of NGC 4258. The two 
serendipitous sources, Q1218H-472 and the putative galaxy clus- 
ter, wiU be addressed in Section 5. 

3.1. Fitting the soft X- my spectrum 

3.1.1. The inclusion ofROSAT data 

M94 showed that the soft X-ray spectrum of NGC 4258 is 
complex, with components arising from thermal plasma emis- 
sion and, possibly, the integrated emission of the X-ray bi- 
nary population. Since we do not consider ASCA data below 
0.8 keV, spectrum models of this region will be poorly con- 
strained with many different models able to explain the soft 
spectrum. To break these degeneracies, we included an archival 
ROSAT PSPC spectra in our analysis. We chose to use the 
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Fig. 1 . — Panel (a) shows the ASCA-GIS2 image. These data have been slightly smoothed for display purposes using a top-hat function with a full-width of 
1 arcmin. Panel (b) shows X-ray contours (taken from the smoothed GIS2 image) overlaid on the optical image from the digitized sky survey. The optical image 
was obtained from the skyview database situated at the NASA Goddard Space Flight Center. 



longest single ROSAT PSPC integration of NGC 4258. This 
observation, which was performed on 1990-Jun-l, has about 
25 ksec of good exposure time. This dataset was obtained from 
the HEASARC archive at the NASA Goddard Space FUght Cen- 
ter and reduced using the FTOOLS routine XSELECT v1.4b. 
ROSAT-PSPC data were used in the range 0.2-2 keV. Good 
agreement was obtained between the ASCA SIS/GIS and the 
ROSAT PSPC in the overlap band between 1-2 keV, thereby al- 
leviating concerns that our analysis will be severely affected by 
poor ROSAT-ASCA cross-calibration. 

3.1.2. Characterizing thermal plasma emission 

Table 1 details our spectral analysis of NGC 4258 (excluding 
our detailed analysis of the iron Une which will be addressed 
in Section 4). As is evident from our spectrum and the previ- 
ous work of M94, the spectrum has both hard and soft com- 
ponents. The simplest model that we attempted to fit consists 
of an absorbed power-law with an additional bremsstrahlung 
component, all absorbed by the Galactic column density of 
^Gai = 1.2 X 10^" cm~^ . This was a dreadful fit to the spec- 
trum, giving x^/dof = 3169/1507. Large residuals exist with 
the model under-predicting the data in the 0.7-1.5 keV range. 
Since this is the energy range in which powerful line emission 
can occur from a thermal plasma with kT ~ 1 keV, we next re- 
placed the bremsstrahlung components with the thermal plasma 
model MEKAL (Mewe, Gronenschild & van den Oord 1985; Ar- 
naud & Rothenflug 1985; Mewe, Lemen & van den Oord 1988; 
Kaastra 1992). Initially, we consider a thermal plasma model 
in which all metals are assumed to have the same fractional ele- 
mental abundance relative to the cosmic abundances of Anders 
and Grevesse (1989). The inclusion of the thermal plasma com- 
ponent (model- A in Table 1) leads to a dramatic improvement 
in the goodness of fit with x^/dof = 1925/1507. However, 
several hne-like residuals in the X-ray spectrum (including one 
at iron Ka energies), as well as a general curvature of the spec- 
trum, prevent this model from being an adequate fit to the data. 
We now discuss extensions of this basic spectral model which 



can adequately explain the observed spectrum. 

3.1.3. An additional Bremsstrahlung-like component 

The curvature in the spectrum requires us to consider an 
additional continuum component. M94 include an additional 
bremsstrahlung component in order to model the diffuse and 
off-nuclear emission in this object. Including a bremsstrahlung 
component (model B) accounts for this curvature and leads to a 
dramatic improvement in the goodness of fit (compare models 
A and B; Ax^ = 298 for 2 additional degrees of freedom). This 
additional component is required even if we allow for abun- 
dance variations in the thermal plasma component (compare 
models B and C; also see below). For completeness, we also in- 
vestigated the possibihty that this additional continuum compo- 
nent is a power-law that does not suffer any intrinsic absorption 
(model D). Scattering of the AGN power-law emission around 
the absorbing material, or non-thermal spatially-extended emis- 
sion from the galaxy or jet would be possible sources of such 
a component. This provides a marginally worse fit than the 
bremsstrahlung-based model, and impUes large scattering frac- 
tions (/scat = 0.16) or powerful distributed non-thermal emis- 
sion (Lx.dist ~ 10^° ergs~^ in the 0.5-lOkeV range). Hence, 
we think this power-law alternative to be unhkely. 

3 . 1 .4. Constraints on the plasma abundance 

Whether this additional continuum is modeled as a 
bremsstrahlung or a power-law component has little influence 
on the best-fit parameters for the direct AGN component (i.e. 
it does not significantly affect the photon index or normaliza- 
tion of the AGN power-law emission, nor the inferred absorb- 
ing column through to the power-law source). However, it does 
effect the best-fit abundances of the soft thermal plasma com- 
ponent (compare models C and B in Table 1). In this sec- 
tion, we attempt to constrain the abundances of the thermal 
plasma under the assumption that this additional component has 
a bremsstrahlung form. We also make the assumption that the 
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Table 1 

Joint ROSAT-FSFC/ASCA spectral fitting results for NGC 4258 



Model 


Components 




r 






ZFe,Ni 


Zlight 




<7 Wpe 


X^/dof 






(xlO^^cm-^) 




(keV) 


(keV) 






(keV) 


(keV) (eV) 






(NH(PO)+. . .) 




















A 


MEKALi 


3.8 ± 0.3 


1.36 ± 0.07 




0.61 ± 0.02 


0.09 ± 0.01 


0.09 ± 0.01 






1930/1506 


B 


BREMS+ 
MEKALi 


9.7 ± 2.0 




> 6.2 


0.36 ± 0.02 


"■■'"-0.06 


„ o„ + 0.08 
"■■^"-0.06 






1627/1504 


C 


MEKALa 


4.3 ± 0.3 


1.44 ± 0.08 




0.56 ± 0.01 


0.10 ± 0.01 


0.35 ± 0.05 






1784/1505 


D 


PO2+ 
MEKAL2 


8 9+°* 

°-^-0.7 


, Qo+0.06 

±.00 Qg 






1.39+_li 


, 1,+S.O 

^•^•^-0.6 






1634/1504 


e 


BREMS+ 
MEKAL2 


9 5+^1 

^•^-0.9 


1.86l°:f3 


> 5.0 


„ ,(,+0.03 
"•^O_0.02 


0.31+_li° 


0.30 ±0.08 






1627/1503 


F 


BREMS+ 


8.2 ±0.9 


^- '^-0.11 


> 4.7 


-0.09 


0.23 ±0.05 


Zc < 2.7 






1557/1503 




MEKAL3 








Zn < 1.6 
Zo < 0.14 
Zno < 0.03 
ZNa < 0.50 

Zms = o.36j;°-Jg 

Zai < 1.3 

zs. = o.7+;,-« 

- J---^_o.8 
ZAr < 1.95 
Zca = 5 ± 2 








G 


BREMS+ 
MFIfAI „-!. 


9 5+^-^ 
•^■^-0.9 


l-°°-0.13 


> 5.0 


U-3D_0.02 




0.30 ± 0.08 


6 45+"-^" 


< 0.20 107+37 


1609/1500 



GAU 



Note. — Model abbreviations follow: NHgai, absorption by the Galactic column density of 1.2 X 10'^" cm^^ ; NH, absorption by a neutral column density of A'H; 
PO, power-law component with photon index F; PO2 second (unabsorbed) power law component with amplitude ratio N1/N2 = 0.16; BREMS, bremsstrahlung 
emission with temperature kT^,.; MEKALi, thermal plasma emission model (Mewe, Gronenschild & van den Oord 1985; Amaud & Rothenflug 1985; Mewe, Lemen 
& van den Oord 1988; Kaastra 1992) with temperature feTpi and all elemental abundances locked together to be a fixed fraction Z of the cosmic value (defined 
by Anders & Grevesse 1989); MEKAL2, thermal plasma emission model in which the relative abundances of Iron and Nickel are fixed together (with the value 
•^Fe.Ni) but can vary independently from the relative abundance of the lighter elements (.Znght); MEKAL3, similar to MEKAL2 except that the individual elemental 
abundances of the lighter elements can vary independently; GAU, Gaussian iron emission line component with rest-frame energy iJune, standard deviation a, and 
equivalent width VFpe- Errors are quoted at the 90 per cent level for 1 interesting parameter (Ax^ = 2.7). 
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a) Spectrum of NGC4258 b) Iron line region 
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Fig. 2. — Panel (a) .shows the joint ROSAT/ASCA spectrum for NGC 4258 fit with model-F from Table 1. The ROSAT daXa are the plain crosses. For clarity, only 
data from the SISO (open squares) and GIS2 (filled circles) instruments on board of ASCA are shown. Panel (b) shows iron line region of the spectrum with less 
severe binning, also referenced to model-F from Table 1 . The presence of a fluorescent iron line is clear. These are both folded spectra (in the sense that they include 
the instrumental response). 



thermal plasma emission can be characterized by a single tem- 
perature. We make these assumptions here in order to be able 
to make progress with these data. Both of these assumptions 
may be invaUd. We must await future high-resolution, high 
signal-to-noise data in order to test and relax these assumptions 
through the use of direct emission line diagnostics. 

We investigated the plasma abundances in a two step pro- 
cess. Firstly, the metals were split into two classes with iron 
and nickel in one class, and all of the Ughter metals in the 
other class. The relative abundances were fixed within each 
class, but the relative abundance of each class was allowed to 
vary independently. This had no effect on the goodness of fit, 
with the best-fit relative abundances of each class being very 
similar (compare models B and E). Since soft X-ray Une-like 
residuals still persist in these fits, the relative abundances of 
all of the light elements were then allowed to vary indepen- 
dently. This leads to a further improvement in the goodness 
of fit (compare models E and F; Ax^ = 70 for 10 additional 
degrees of freedom). After fitting model F (which includes an 
additional bremsstrahlung component; see below) it can be seen 
from Table 1 that the abundances of C, Ne, Na, Al, and Ar are 
poorly constrained with only weak upper Umits possible. This 
is due to the lack of strong emission lines from these elements 
in the well-calibrated region of the A5CA-S1S/G1S. On the other 
hand, the abundances of Mg and Fe are well constrained (with 
Zug = 0.36 and Zpe = 0.23) due to the detection of fairly 
strong emission fine complexes associated with these elements. 
It is worth noting the apparent extreme overabundance of Cal- 
cium (Zcii = 5 ± 2). This result must be viewed with sus- 
picion since the dominant Calcium emission lines emerge just 
below our usable ASCA band (which has a lower-energy cut- 
off at 0.8 keV for the CIS and 1 .0 keV for the SIS). Hence, any 
sUght ASCA/ROSAT cross-caUbration problem might be mani- 
fested as an extreme Calcium abundance. 

3.1.5. Fluxes and luminosities 



Using model F from Table 1, and assuming a distance of 
7.2 Mpc, we can compute the absorption-corrected luminosi- 
ties in each spectral component. The result is shown in Table 2. 
It can be seen that the power-law component dominates the en- 
ergetics of the X-ray band by almost an order of magnitude. On 
energetic grounds, the thermal plasma emission can be pow- 
ered via the absorption and re-emission of 10-25% of the hard 
power-law component (depending upon how low in energy the 
power-law component extends). Our data cannot probe the na- 
ture and origin of this thermal emission beyond pointing out 
the basic energetics. However, future observations with Chan- 
dra and XMM will be able to study the spatial distribution of the 
various spectral components that we have noted. One will then 
be able to address whether the thermal emission is located in 
the immediate vicinity of the AGN (as, for example, if it orig- 
inated from an accretion disk wind) or distributed on scales of 
100 pc or greater (as in the case of a galactic superwind). 

3.2. Characterizing the iron Ka emission line 

Residuals due to the iron fine are clearly visible when mod- 
els A-E are fit to these data (see Fig. 2b). We choose to use 
model D of Table 1 as a base model for our iron line inves- 
tigation. We do not use model E (i.e. the best fitting model 
from Table 1) because the excessive number of free parameters 
makes the iron line error analysis unnecessarily difficult. We 
have verified that our iron line results are not affected in any 
significant manner by the choice of using model D as a base 
model rather than model E. We initially fit this line by adding a 
Gaussian emission feature (model F of Table 1). The improve- 
ment in the goodness of fit is significant at more than the 90 per 
cent level (Ax^ = 18 for 3 additional degrees of freedom). The 
full-width half maximum (FWHM) of the line is constrained to 
be less than 22000 km , and is consistent with being zero 
(i.e. the line is unresolved). The inferred centroid energy of 
E = 6.45lo;o7 consistent with the Ka Une energy from 

iron with an ionization state of less than Fexvii. The equivalent 
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Table 2 

x-ray fluxes and luminosities for ngc 4258 based on model f of table 1. the bolometric luminosity of the thermal 
plasma emission can been calculated under the assumption that a single temperature describes this material. all 
luminosities assume a distance of 7.2 mpc and have been absorption-corrected (i.e. both galactic and intrinsic 

absorption have been removed). 



parameter value 

observed 0.5-10 keV flux 2.1 x 10~^^ ergcm"^ 

observed 2-lOkeV flux 5.8 x 10^^^ crgcm"^ s~-^ 

power-law 0.5-10 keV luminosity 7.9 x IC^^ergs"^ 

thermal plasma 0.5-lOkeV luminosity 1.1 x 10*°ergs~^ 

thermal plasma bolometricluminosity 2.1 x lO^^ergs"^ 

bremsstrahlung 0.5-10 keV luminosity 9.2 x 10"^^ergs~^ 



width of the Hne is Wfc = lOT+ly cV. 

We discuss the origins of this emission hne and the possible 
impUcations for the nature of the accretion disk in Section 4. 

3.3. The long term variability of NGC 4258 

ASCA observed NGC 4258 on 4 previous occasions, thereby 
allowing us to examine the variability of this AGN on 
timescales of several years. We retrieved all of the available 
ASCA data on NGC 4258 from the HEASARC public database 
located at the NASA Goddard Space Fhght Center. These data 
were reduced in the same manner as described in Section 2, 
and then fitted with spectral model B from Table 1 . Since these 
archival data possess significantly lower signal to noise, the 
temperatures of the thermal plasma and bremsstrahlung com- 
ponents were fixed to the best fit values from Table 1. The 
normalization of the bremsstrahlung component was also fixed. 
Thus, this investigation is probing changes in the AGN power- 
law and absorbing column density. Note that only GIS data 
were included for the 5-May-93 observation, since NGC 4258 
faUs very close to the SIS chip boundaries in this observation. 

Table 3 summarized the results of this study. The photon 
index of the power-law emission is consistent with being con- 
stant over the 6 years covered by these data. However, there 
are indications of variations in both the absorbing column den- 
sity (which seems to decrease by ~ 30% between the 1993 
and 1996 observations) and the 5-lOkeV flux (which almost 
doubles from the 1993 to the 1996 observations, and then drops 
back to the 1993 level by the time of the 1999 observation). The 
flux in the 5-10 keV range is dominated by the AGN power-law 
component and is little affected by the absorbing colunm. 

4. IMPLICATIONS OF THE IRON LINE IN NGC 4258 

4. 1 . Pure accretion disk models for the iron line 

The interest in the iron hne lies in its ability to diagnose 
the nature of the accretion disk. Suppose that the only sig- 
nificant source of iron Ka line emission in NGC 4258 is the 
AGN accretion disk. We can then model the observed iron 
hne as being from the surface of the accretion disk around a 
non-rotating (Schwarzschild) black hole by using the DISKLINE 
model within the XSPEC software package. We will make the 
assumption that the disk is flat inside of the masing radii and so 
set the inclination of the inner disk to be i — 85°. The outer ra- 
dius of the iron line emitting region is set to rout — 10^ GM/c? 
(the radius of the maser disk; note that the iron hne fits are very 
insensitive to the actual value of rout provided it is sufficiently 



large). The line energy was fixed at the value appropriate for 
Ka emission from weakly ionized iron, E — 6.40 kcV. Free 
parameters in the fit are the inner radius of the line emitting 
region rbr, the index describing hne emission as a function of 
radius /? (where the surface emissivity e oc r^), and the normal- 
ization of the line. Figure 3a reports the confidence contours 
that result when this accretion disk model is fit to the iron line. 
It can be seen that the narrowness of the line requires either flat 
emissivity as a function of radius (/? > —2) or an inner edge to 
the line emitting region at greater than a hundred gravitational 
radii. 

If the X-ray emission traces the viscous dissipation in the 
disk, and the disk is geometrically-thin and radiatively-efficient 
so that it can be described by a Novikov & Thorne (1973) 
model, the emissivity index should tend to /3 = —3 outside 
of the inner 20 gravitational radii or so. In this case, and given 
the assumption that the observed iron line is from the accretion 
disk, we are led to the conclusion (at the 90% confidence level) 
that the line emitting region has an inner edge at ^ 100 GM jc? . 
Such an inner edge may correspond to the point at which the 
disk surface becomes ionized, or where the disk undergoes a 
transition to a hot (possibly advection dominated) state. On the 
other hand, if the X-ray source is in the form of a geometrically- 
thick corona with size D, the emissivity index will be fairly flat 
for r < D, and will tend to (3 = — 3 for r >> D. In this case 
(and again, with the assumption that the observed iron line is 
from the accretion disk), we must conclude that the corona is 
large D ^ WOGM/c^. 

4.2. Hybrid disk/non-disk iron line models 

An alternative that we must explore is one in which a some 
fraction of the observed narrow iron line is produced by dis- 
tant material not directly related to the accretion disk. In this 
case, any broad emission line from the accretion disk would 
be blended with this narrow hne and, possibly, buried in the 
noisy continuum spectrum. To investigate this possibility, we 
suppose that a narrow iron line from some non-disk origin con- 
tributes to the observed spectrum with an equivalent width of 
WKa,nar- Fixing the line emissivity profile of the disk to the 
canonical (3 = —3 case. Fig. 3b shows the confidence con- 
tours on the WKa.nar — plane. It can be seen that the line 
emission needs an inner edge or break at rbr ^ 50GM/c^ un- 
less the additional narrow hne source contributes at the level 
;^40eV. 

We cannot rule out, in any rigorous sense, the possibility that 
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Table 3 

Historical A5CA data for NGC 4258. Spectral parameters shown here are derived by fitted model B of Table 1 to 
THESE data. See text for further details. Errors are quoted at the 90 per cent level for 1 interesting parameter 

(Ax' = 2.7). 



Obs. 5-May-93 23-May-96 5-June-96 18-Dec-96 15-May-99 



MJD 


49122 


50226 


50239 


50435 


51313 


SIS Exposure (1000 s) 


39 


23 


31 


26 


170 


r 


1 7O+0.22 


1 7-1+0.18 

^- < '^-o.n 


1.83 ±0.13 


1.87 ±0.15 


1 OO+0.404 


A^hCIO'^ cm-2) 


13.6tli 


9.2 ±0.9 


°-°-0.6 


9.7 ±0.8 


Q ^+2.1 


GIS2 5-lOkeV flux (IQ-^^ ergcm-^ s'^ ) 


5.1 


8.3 


8.8 


9.5 


4.0 


X^/dof 


652/511 


639/542 


909/728 


828/632 


1627/1504 



a) ^ b) 




Fig. 3. — Panel (a) shows the confidence contours on the /3-rbr plane when the observed iron line is modelled as arising from an accretion disk in which the line 
emissivity is e = for r < r]-,r, and e oc for r > rt^. Panel (b) shows the confidence contours on the WKa,nar-^bY plane, when it is assumed that there is 
also a narrow iron line from some non-disk source with an equivalent width WfCa nar- The line emissivity of the disk is fixed to the canonical values of /3 = — 3 in 
this panel. In both figures, the accretion disk is assumed to have an inclination of i = 85°. The confidence contours are at the 68% (dashed line), 90% (soUd line) 
and 95% levels (dotted line) for two interesting parameters. In both plots, the lower-left is the forbidden region and the upper-right is the allowed region. 



most (or all) of the observed iron line originates from matter 
that is not associated with the accretion disk. However, sim- 
ple arguments lead us to disfavor such a scenario. Consider 
iron Une emission in a geometrically-thick torus surrounding 
the accretion disk of NGC 4258. An upper limit to the column 
density of this structure along our line of sight to the AGN is 
given by the observed column density of A^h ~ 10^'^ cm~^ . If 
we suppose that this torus is in the same plane as the accretion 
disk (so that we are also viewing it edge-on), it is plausible to 
assume that we are looking through the optically-thickest part 
of the torus. By considering the case in which the torus com- 
pletely surrounds the AGN with uniform coluimi density along 
aU radii, an upper Umit to the equivalent width of the iron line 
is given by: 

WFe.max = .E2^e^FeA^H^abs,Fe ^ dE , (1) 

where E'une = 6.4 keV is the energy of the emission line, 
Fpe = 0.33 is the fluorescent yield of the transition, Zabs.Pe ~ 



4 X 10"'' is the density of iron relative to hydrogen, and a{E) is 
the energy dependent photoionization cross section for ioniza- 
tion from the Is shell. Here, we have approximated the photon 
index of the ionizing AGN power-law to F = 2 and assumed 
that the central X-ray source is isotropic. Using the photoion- 
ization cross-sections for neutral iron of Vemer & Yakovlev 
(1994), this yields 

.max 

«65eV. (2) 

However, there are several reasons why this upper limit 
would almost certainly not be achieved. Firstly, modeling of 
the accretion disk warp strongly suggests that our line of sight 
intercepts the disk and that the bulk of the colunm density which 
obscures the AGN originates in the disk (Hermstein, priv. com- 
munication). Therefore, we might expect significantly smaller 
column densities along lines of sight that have smaller inclina- 
tions angles relative to the accretion disk. Secondly, the accre- 
tion disk may well occult half of this fluorescing cloud, thereby 
reducing this prediction further. Thus, the true iron line from 
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surrounding non-disk material may well be reduced from our 
naive prediction by a factor of several. 

A high column density torus that is misaligned with the ac- 
cretion disk so as to leave the X-ray unobscured is still a viable 
source for the observed narrow iron line. Such a torus must ex- 
ist on size scales significantly larger than the masing accretion 
disk, or else it would hydrodynamically disturb the disk and 
lead to strong deviations from the observed Keplerian rotation. 

4.3. Equivalent width limits on a "Seyfert-like" broad iron 

line 

Many higher luminosity Seyfert galaxies display iron Unes 
which are so broad that the line emitting region is thought to 
extend down to near the radius of marginal stability (dGM/c^ 
for a non-rotating black hole). The typical emissivity index (3 
Ues between —2 and —3 (Nandraet al. 1997). As shown above, 
if WKa,nar ^ 40 eV then no such relativistically broad com- 
ponent is required by our fits to NGC 4258. However, even 
when all of the observed line is assumed to arise from some 
other structure, we can still obtain an upper Umit to the equiv- 
alent width of any such "Seyfert-like" component by fixing 
rbr = QGM/c^ in the spectral fitting, and including an ex- 
plicit narrow Gaussian component to fit the observed line. The 
resulting upper limit on the equivalent width of the relativistic 
iron line, as a function of the assumed line emissivity index, is 
shown in Fig. 4 (thick solid line). 

If NGC 4258 possesses an iron line emissivity inner accre- 
tion disk similar to higher luminosity Seyfert galaxies, one 
might naively expect that hmb-darkening due to absorption in 
the outer layers of the accretion disk will reduce the equiva- 
lent width of any such broad iron line to very small values (e.g. 
George & Fabian 1991). From the work of Ghisellini, Haardt 
& Matt (1994), the hmb-darkening in the plane-parallel case is 
well described by the expression 



WKaje = 0) 

In 2 



COS ^ In I 1 



1 



cos 6* 



(3) 



where 9 is the inchnation angle. Nandra et al. (1997) shows that 
most Seyfert 1 galaxies have a broad iron line with Wxa in the 
range 200-400 eV. Taking these to be representative of the face- 
on values, an inchnation of 9 = 85° and this limb-darkening 
law reduces the expected equivalent width to 66-1 32 eV. This 
is below our detection threshold. 

Two important effects modify this estimate. Firstly, relativis- 
tic light bending means that significant portions of the inner- 
most regions of the disk are viewed significantly more face- 
on than otherwise might be thought. The dotted and dashed 
lines on Fig. 4 show the effects of light bending on the pre- 
dicted broad line equivalent width (assuming a face-on equiva- 
lent width, WKa{9 = 0), of 200 eV and 400 eV respectively). 
These lines have been computed using the relativistic code pre- 
sented in Reynolds et al. (1999) combined with the above hmb- 
darkening law. Figure 4 shows that the broad line would not be 
detectable unless WkJ9 = 0) > 300 cV and /3 > -3. 

Secondly, the above limb-darkening law assumes that the 
T = 1 surface of the accretion disk is strictly planar. It 
is very unlikely that this will be the case. MHD turbulence 
and violent plasma processes within the disk corona will in- 
evitably corrugate this surface thereby reducing the effects of 

'Note that, following the convention used in xspec, we have defined the equivalent width with respect to the continuum level at the energy where the line peaks in 
photon flux. For centrally concentrated disk illumination (i.e. very negative the iron line peaks at 8 keV. Since the continuum flux is a strongly declining function 
of energy, the equivalent width of such a line can exceed the face-on value. 



limb-darkening for highly inclined sources. In the extreme 
case, most of the matter in the outer layers of the accretion 
disk may be concentrated into dense filaments or clumps (e.g. 
see the MHD simulations for the 'Channel solutions' case pre- 
sented by Miller & Hawley 2000), thereby removing orienta- 
tion and limb-darkening effects almost entirely. The dot-dashed 
and the dot-dot-dot-dashed lines on Fig. 4 show the predicted 
equivalent width in the case of no hmb-darkening, assuming 
WKa{9 = 0) = 200 eV and 400 eV, respectively^. In this lim- 
iting case, we would expect to be sensitive to such a broad line 
for all Seyfert-hke values of /3 and WKa{9 = 0). 

To summarize, these data just begin to constrain the inter- 
esting regions of parameter space for any relativistic iron line 
component in NGC 4258. There are two possible reasons for 
our non-detection of a broad iron line. Firstly, such a line might 
be genuinely absent thereby setting NGC 4258 aside from its 
higher luminosity counterparts. Secondly, if NGC 4258 pos- 
sesses a relativistically broad iron fine with WKa{6 = 0) ^ 
300 eV and limb-darkening is important, then we would not 
expect to detect this line with these data. However, we can 
rule out the presence of a relativistically broad iron fine with 
Wkcx{0 = 0) ^ 300 eV and /? > -3. Neither do these data al- 
low a Seyfert-like line in which limb-darkening is unimportant. 

5. THE SERENDIPITOUS SOURCES 

5.1. Q1218+472 

The point source 7arcmins to the west of NGC 4258 is 
readily identified with the = 0.40 quasar Q1218H-472. 
By choosing a slightly smaller extraction radius than normal 
(2.5 arcmins), we can isolate the GIS counts from this source 
with relatively little contamination from NGC 4258. The re- 
sulting GIS background subtracted count rate is 4 x 10~^ cps 
per detector. No temporal variability was observed in the re- 
sulting light curve, although the low count rate results in weak 
limits on possible variability. This source is not within the field 
of view of the SIS. 

The two GIS spectra were fitted with a model consisting of a 
power-law subjected to both Galactic absorption (with a col- 
umn of -/Vcai = 1.2 X 10^°cm~^) and intrinsic absorption 
with the redshift of the quasar. The resulting best fit parameters 
are T = lAQtoW and iVn < 5 x lO^i cm-^ (see Fig. 5 for 
the confidence contours that result from this fit). The (observer 
frame) 2-lOkeV flux is Fz^io = 4.0 x IQ-^^ ej-gcm^^ g-i . 
Assuming a redshift of z = 0.40, a Hubble parameter of 
Hq = 65kms~^ Mpc^^, and an acceleration parameter of 
go = 0.5, the rest-frame 2-10 keV luminosity of this source is 
i'2-io = 1.6 X 10^^ erg s~^ . In computing the flux and lumi- 
nosity, we have applied a correction factor of 1.5 to account for 
the smaller than usual extraction radius. This correction factor 
was estimated by extracting counts from other ASCA observa- 
tions of bright Seyfert galaxies with various extraction radii. 

While the flux and luminosity estimate are robust, the spec- 
tral results should be treated with caution since the GIS count 
rates for this source are below the threshold at which the ASCA 
Guest Observer FaciUty recommends spectroscopy. Combin- 
ing our X-ray spectrum with the optical spectrum of Burbidge 
(1995), the spectral index between the (rest-frame) 2500A and 
2keV emissions is ctox ~ 10. Thus, Q1218-H472is relatively 
X-ray bright given its optical flux. 
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Fig. 4. — Constraints on the presence of a "Seyfert-like" iron line in the case where the observed narrow line is modelled by a separate narrow Gaussian compo- 
nent. The broad line is modelled as originating in an accretion disk around a Schwarzschild black hole, with a disk inclination of i = 85° and an inner line emitting 
radius of r^r = 6GM/c^. The solid line shows the upper limit on the equivalent widths a function of the emissivity index 8. The dotted and dashed lines show 
the theoretical expectation, taking into account limb-darkening and light bending effects, assuming that the iron line has an equivalent width of 200 cV and 400 cV, 
respectively, when the disk is viewed face on. The dot-dashed and dot-dot-dot-dashed lines show the theoretical expectation when limb-darkening is absent (see 
text) assuming that the iron line has an equivalent width of 200 eV and 400 eV, respectively, when the disk is viewed face on. 
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Fig. 5. — Constraints on the intrinsic column density and photon index for Q1218-H472. Shown here are the la, 90 per cent and 95 per cent confidence contours 
for two interesting parameters. 



5.2. The north-west source: a putative galaxy cluster 

The northwest source aligns well with the extended ROSAT- 
PSPC source found by Pietsch et al. (1994; hereafter P94). P94 
note a concentration of galaxies on a deep optical plate. By not- 
ing that the brightest galaxy is approximately 19 magnitude in 
the visual band, they estimate the redshift of the cluster to be 

Using standard extraction radii, we have extracted the GIS 
spectra and light curves for this source (note that this source. 



too, fell outside the SIS field of view). The GIS count rate was 
7 X 10""^ cps per detector. Given the possible identification of 
this source with a cluster of galaxies, we fit the spectrum with 
a thermal plasma model (MEKALi of Table 1) modified by 
Galactic absorption. We also allow the redshift of the cluster to 
be a free parameter. The best fit plasma temperature, relative 
abundance, and redshift are kT = S.ltJJ keV, Z = 0.49l^;32 
and z = 0.28 ± 0.05 respectively. Figure 6 shows the two GIS 
spectra along with this best fit model. It can be seen that the 
model is a good fit to these data (with x^/dof=202/217). This 
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a) Putative galaxy cluster 
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Fig. 6. — GIS spectrum of the north-west object, probably a 2 ~ 0.3 cluster of galaxies. The best fit thermal plasma model is shown (see text for details). 



X-ray determination of the cluster redshift relies on fitting the 
cluster iron emission line to a bump in the GIS spectrum at 
^ 5 keV. Since instrumental GIS features are possible at these 
energies, this must be taken as tentative. 

The observer-frame 2-lOkeV flux of this source is -F2-10 = 
6 X 10^^^ ergcm^^ s^^ . Assuming a cluster redshift of z = 
0.28, a Hubble parameter of Hq = 65kms~^ Mpc~^, and an 
acceleration parameter of qo = 0.5, the rest-frame 2-10 keV 
luminosity of this source is Z/2-io = 1-4 x lO^^ergs"^. 
These results are completely consistent with the known re- 
lationship between cluster temperature and X-ray luminosity 
(e.g. Markevitch 1998). If, instead, this source is at the 
same distance as NGC 4258, the 2-10 keV luminosity is only 
L2-in = 4x 1039- 



erg s 

6. DISCUSSION AND CONCLUSIONS 



The X-ray spectrum of NGC 4258 shows 4 distinct compo- 
nents: 

1. a power-law component with a 0.5-lOkeV luminosity of 
7.9 X lO'^Oergs-i and photon index T = UgtiH- 
This component is readily identified as the primary X- 
ray emission from the AGN central engine itself. This 
emission is absorbed by an intrinsic column density of 
iVn = 8.2 ± 0.9 X 10^2 cm"^ ^ ^ote that this absorb- 
ing column is a factor of 1.4 smaller than that obtained 
by M94. Also note that our value of Nu is insensitive to 
exactly which of the acceptable spectral models we use 
(cf. Gammie, Narayan & Blandford 1999; Chary et al. 
2000). 



thermal plasma emission from optically-thin gas at a tem- 
perature of kT ~ 0.5 keV. The total luminosity of this 
component is ^ 2 x IC^^crgs^^, approximately 10- 
25% of the power-law luminosity (depending on how low 
in energy the power-law extends). 



3. a bremstrahlung component with kT > 5keV. This 
component may represent thermal emission from a very 



hot gaseous component in the anomalous arms and/or the 
integrated hard X-ray emission of the X-ray binary pop- 
ulation. 

4. a narrow fluroescent Ka emission line of cold iron with 
an equivalent width of Wku = 1071 gy eV. 

We suggest that the bulk of the narrow iron emission line 
originates from the accretion disk. With this assumption, we 
have shown that the most of the line emission must originate 
at relatively large radii (r ^ lOOGM/c^) in order to produce 
the small line width. This, in turn, implies a large X-ray source 
(with a size D lOOGM/c^), in contrast to the situation nor- 
mally found in higher-luminosity Seyfert galaxies (e.g. Nandra 
et al. 1997). Note that a model in which the X-ray source is 
small (D ~ lOGM/c^) but the inner region of the disk are too 
hot/ionized to produce line emission has difficulty producing 
the observed equivalent width of the line. Thus, under the as- 
sumption that the bulk of the observed iron line originates from 
the accretion disk, there appears to be a clear difference in the 
size/structure of the X-ray source between this low-luminosity 
source and higher-luminosity Seyfert galaxies. 

However, we cannot rule out the possibility that the narrow 
iron line originates from some previously undetected distant 
matter not associated with the accretion disk. If the observed 
line does originate from a non-disk source, our data are con- 
sistent with but do not require the presence of a "Seyfert-like" 
relativistic broad iron Une. Even in this case, interesting con- 
straints can be placed on the parameter space of possible rela- 
tivistic broad iron lines. 

This AGN is ripe for study with the new generation of X- 
ray observatories. The high spatial resolution of Chandra will 
allow us to study large scale structure in this source directly. 
We will be able to pinpoint the location and nature of the ther- 
mal component seen in our ASCA data. We will also be able 
to spatially separate emission from the anomalous arms, the X- 
ray binary population, and the AGN. High signal-to-noise spec- 
troscopy with XMM (and, in the longer term, Constellation-X) 
will allow the iimer accretion disk to be probed. 
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